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Toroidal carbon nanotubes can serve as hosts for encapsulated loops of atomic metal wires. Such 
composite structures have been analyzed using density functional theory for a semiconducting C120 
torus encapsulating chains of Fe, Au and Cu atoms. The sheathed metal necklaces form a zigzag 
structure and drops the HOMO/LUMO bandgap to less than 0.1 eV. The iron composite is ferro- 
magnetic with a magnetic moment essentially the same as that of bcc iron. The azimuthal symmetry 
of these toroidal composites suggests that they may offer novel elecromagnetic properties not asso- 
ciated with straight, metal-encapsulated carbon nanotubes. 



INTRODUCTION 

Fullerene-based toroids combine the novel characteris- 
tics of carbon nanotubes with electromagnetic and me- 
chanical properties unique to closed-ring geometries. An 
intriguing possibility, examined in this paper and shown 
in Figure [l] is that atomic chains of metal can be encap- 
sulated within such toroidal carbon nanotubes (TCNs). 
The carbon torus would serve as a stabilizing sheath for 
the atomic metal necklaces resulting in rigid, stable struc- 
tures at finite temperatures with novel paramagnetic or 
ferromagnetic properties properties. Beyond applications 
that exploit the properties of independent rings, assem- 
blies of such ferromagnetic tori might be expected to self- 
assemble into columnar structures with unusual macro- 
scopic behavior. The electronic structure and mechanical 
stability of these robust atomic necklaces are investigated 
using density functional theory (DFT) with the objective 
of establishing ground state geometries, binding energies 
and magnetic moments. 

At least four types of single-wall fullerene toroids have 
been considered since the idea was first proposed fifteen 
years ago[l . Two variations of Type 1 TCNs exploit 
the curvature induced by pentagonal and septagonal de- 
fects to generate polygonal graphitic tori pQ, [2] and rings 
as small as C120 have been considered [3]. The sym- 
metric Type 2 TCNs are toroidal poly hexes [4] created 
by joining the ends of single- wall carbon nanotubes [5]. 
Perhaps the most technologically promising category of 
toroid is derived from Cqq molecules that are joined with 
pentagonal and septagonal rings and looped into undu- 
lated rings [6 ,[7 . These are referred to here as Type 
3 toroids. Finally, Haeckelite tubules [8 , based on a 
two-dimensional bravais lattice with 5, 6 and 7- member 
rings, can also be joined at the ends to make the Type 4 
TCNs [6], 0. The first two types of TCNs are theoretical 
constructs that are thought to be mechanically stable. It 
has been estimated that TCNs of Type 1 are in fact more 
stable than Ceo molecules [3 , while Type 2 TCNs are ex- 
pected to be stable provided they have a sufficiently large 
toroidal diameter [5], 0. Carbon toroids of Types 3 and 
4, on the other hand, have been realized experimentally 
and are demonstrably st able [6], [7]. 
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FIG. 1: Atomic chains of Fe atoms encapsulated within three 
types of toroidal carbon nanotubes. A bucky ball is shown at 
lower right for size comparison. All structures are drawn to 
scale. 



For all toroids, delocalized tt electrons can be in- 
duced into ballistic, persistent motion by the applica- 
tion of a magnetic field perpendicular to the plane of the 
torus. The resulting orbital magnetic response of tori 
can be substantial and has been investigated using tight 
binding theory with the magnetic vector potential ac- 
counted for using London theory [T0],[TT]. Type 1 toroids 
are predicted to exhibit a large diamagnetic response 
Q2] , Q2] , [5] , while a class of conducting Type 2 TCNs are 
expected to be able to generate colossal paramagnetic 
moments [14] , [5] , [15] . More recently, carbon toroids of 
Types 3 and 4 TCNs have been shown experimentally to 
have paramagnetic behavior consistent with tight bind- 
ing model predictions [6 . The septagonal and pentagonal 
rings in these structures are predicted to result in ferro- 
magnetic behavior at sufficiently low temperatures in the 
absence of any flaws. All types of TCNs are expected to 
exhibit Aharanov-Bohm effects [13 . 

Paralleling activity associated with toroidal fullerene 



2 



structures have been advances in the fabrication and 
analysis of isolated atomic chains of metal atoms. Metal 
contacts created with a scanning tunneling microscope 
can be slowly pulled apart to extrude an atomic chain 
of Au atoms [16], [T7]. Stable monatomic chains of Au 
atoms can support currents of 80 /iA-clear evidence 
of ballistic transport in this alternate manifestation of 
one-dimensional conductors [TS] -and the quantum con- 
ductance in these metallic atomic chains has also been 
measured [1 7 . As subsequent DFT analysis of atomic 
chains of Au atoms, intended to address observed anoma- 
lies, indicated that zigzag structures are more stable than 
a linear chains as a result of reduced transverse kinetic en- 
ergy of electrons [18]. Such configurations are most likely 
metastable since previous investigations have identified 
polyhedral ground state structures for small clusters of 
Fe and Au atoms [19], [20]. 

Iron structures sheathed within straight single wall car- 
bon nanotubes have been studied both computationally 
and experimentally [2T] . [22] . A DFT analysis of atomic 
chains of Fe within single-wall carbon nanotubes con- 
cluded that the spacing between Fe atoms and their 
nearest C neighbors determined the character of the Fe 
wire[2T]. The Fe-C bond length of 1.9 A associated with 
(3, 3) tubes resulted in a composite with semi-conductor 
character and no magnetic moment. On the other hand, 
an Fe-C bond length of 2.1 A associated with (4, 4) tubes 
resulted in a composite with the metallic character and 
magnetic moment of 2.6/i# per atom[2T]. The geometry 
of such tubes precludes the generation of axial persis- 
tent currents and a colossal paramagnetic response, but 
they might be used to help predict the conduction and 
ferromagnetic nature wires encapsulated within TCNs. 

The present work brings together the unusual behavior 
of toroidal carbon nanotubes and atomic metallic chains 
within a simple setting: a Type 1 C120 torus with D 5 ^ 
symmetry [2] encapsulating closed, monatomic chains of 
Au, Fe, and Cu atoms. The iron necklace is of partic- 
ular importance because of the possibility that it has a 
ferromagnetic nature at finite temperatures. The gold 
necklace was also considered, though, because of the fo- 
cus that Au has received in atomic metal wires [T6].[T7]. 
Analysis was Cu atoms was thought to offer an interme- 
diate composite since it is a noble metal, like Au, and a 
Period 4 metal like Fe. The Type 1 TCN chosen for anal- 
ysis is a polygonal C120 torus with D$d symmetry-chosen 
for its predicted stability and minimal size [2 . The C120 
torus is composed of an inner surface of ten septagons and 
an outer surface of ten pentagons patched with twenty 
hexagons as shown in Figure [2] Encapsulated atomic 
loops composed of ten metal atoms were considered in 
order to exploit maintain the computational expedient of 
D^d symmetry while creating M-M bond lengths on the 
order those calculated for isolated metal chains. Den- 
sity functional theory (DFT) was used to determine the 
ground state configurations and to quantify M-M bond 




FIG. 2: Break-out view of the inner and outer structures of 
the Ci2oFeio molecule with shading used to highlight polyg- 
onal defect types. 



length, minimum C-M bond length, magnetic moment, 
density of states, conduction character, binding energy 
of the encapsulated metal atoms. 



METHOD 

All calculations were performed with the real-space, 
numerical atomic orbital, DFT code, DMOL 3 [23 . A 
norm conserving, spin unrestricted, semi-core pseudopo- 
tential approach was employed with electron correlation 
accounted for using the Perdew-Wang generalized gra- 
dient approximation (GGA)[24. Double numerical basis 
sets were used with 4p polarization functions for Fe atoms 
and 3d polarization functions for C atoms. A smearing 
value of 0.272 eV was used for all structures. 

As a check on the method, an Fe dimer was determined 
to have a ( 7 A n ) collinear ferromagnetic ground state with 
a bond length of 2.00 A and magnetic moment of 7.0/xb 
by Hirshfeld analysis-consistent with both experimental 
and other DFT results [25 . As a further check, the atomic 
spacing, binding energy and magnetic moment of bcc Fe 
was evaluated to be 2.48 A, —4.89 eV and and 2.22/i^- 




FIG. 3: Ground state structure for a periodic chain of Fe 
atoms. 



reasonable matches to the experimental values of 2.47 
i, -4.32 eV and 2.22/i B per atom [26], [27], [28]. DFT 
values of 2.76 eV and 2.08/i# per atom were obtained 
elsewhere [29 . 



RESULTS 
Free Standing Metal Chains 

Periodic structures of pairs of free standing Fe, Au 
and Cu atoms were analyzed in order to understand the 
character of metallic chains prior to encapsulating them 
in carbon toroids. The ground state configuration for the 
Au chain was found to be a zigzag structure with a bond 
length of 2.83 A and an angle between adjacent atoms of 
57°. There is also a local equilibrium of higher energy at 
134°, and this is close to that obtained by Sanchez-Portal 
et al. who identified an equilibrium angle of 131° [18]. 
Analogous results, summarized in Tables 2 and 3, were 
obtained for Fe and Cu chains, and all three chains ex- 
hibit ground state bending angles of approximately 60°; 
an affinity for symmetric, three-atom bonding is optimum 
within the geometric constraints imposed. As shown in 
Figure [3] such structures amount to two-dimensional, tri- 
angular lattices structures rather than linear chains. Key 
property data is summarized in the tables. All positive 
eigenvalues of the Hessian matrices associated with peri- 
odic cells of four atoms indicate at least linearized stabil- 
ity at K for all three chains. A Jahn- Teller distortion, 
manifested in both bond angle and length, was noted 
for metallic strands of eight atoms or more, and this is 
consistent with previous ab initio investigations [30 . 

Atomic Metal Wires within SWCNT 

Single-strand atomic wires of metal atoms encapsu- 
lated within single-wall carbon nanotubes (SWCNT) are 
of interest on their own and can be used to anticipate 
the electronic structure and stability of wires encapsu- 
lated within carbon tori. The tube diameter of the C120 
torus is between that of (3, 3) and (4, 4) SWCNT, and 
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(3,3)-Fe SWCNT 




FIG. 4: Ground state structures of (3,3) and (4,4) SWCNT 
embedded with atomic Fe wire. 



so these fullerenes were chosen for more detailed anal- 
ysis with metal wires. Previous DFT studies considered 
atomic wires of Fe within SWCNT [21], [31]. The analysis 
for (n, n) tubes considered only one metal atom per unit 
cell, though, and therefore did not allow for the possi- 
bility of the metal chain forming a zigzag configuration. 
Since the toroidal structure should allow for such corru- 
gation, the SWCNT analysis was performed for double- 
cell, periodic structures with both Fe, Au and Cu wires. 
The resulting equilibrium configurations for are shown 
in Figures [4j [5] and [5] with property data summarized in 
Tables 1 through 3. As expected, the encapsulated wires 
do exhibit a zigzag structure although the carbon sheath 
flattens the wires relative to their free standing configu- 
ration. For the (4, 4) tube, the Fe wire causes a diametral 
tube strain of 0.94% while the Au wire causes a larger di- 
ametral expansion with 2.44%. Reassuringly, the Cu wire 
causes a dilation between these two-2.1%. In the case of 
the smaller (3, 3) tube, the diametral expansion of the 
Au tube is 8.7%, while the Fe tube shows a non-uniform 
distortion with bi-lateral strains of 8.2% and —5.2%. The 
Cu tube is uniformly dilated, like Au case, with a strain 
of 3.4%. The allowance for zigzagging within the tubes 
significantly impacts the electronic structure of the com- 
posite. If a periodic cell half as long is considered, so 
that there is only one Fe atom per period, the bandgap 
is 0.87 eV instead of 0.48 eV. The first number is con- 
sistent with the 0.88 eV obtained previously [51]. The 
magnetic moment of both straight and zigzagging chains 
of Fe atoms were extinguished within (4, 4) CNTs. An 
earlier analysis on straight, sheathed chains reported a 
magnetic moment of 2.6/i# per Fe atom[5T]. 
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(3,3)-Au SWCNT 




FIG. 5: Ground state structures of (3,3) and (4,4) SWCNT 
embedded with atomic Au wire. 



(3,3)-Cu SWCNT 




FIG. 6: Ground state structures of (3,3) and (4,4) SWCNT 
embedded with atomic Cu wire. 



C120 Torus 

The C120 torus was first analyzed in the absence of 
atomic wire, and the energy per C atom was found to 
be less than that of Ceo by 0.09 eV (smearing) and 
0.16 eV (no smearing). This is in contrast to an ear- 
lier tight binding calculation which predicted that the 
toroidal state would be of lower energy [2] . The differ- 
ence is small, though, and the C120 torus should be very 
stable. A reactive molecular dynamics code [9] was used 
to estimate that the structure is stable to at least 2000 K. 
The torus has no intrinsic magnetic moment and exhibits 
a bandgap of 0.91 eV between its highest occupied and 
lowest unoccupied molecular orbitals (HOMO/LUMO). 
For comparison, the HOMO/LUMO bandgap of Ceo was 
found to be 1.5 eV-close to the experimental measure- 
ment of 1.6 eV[32]. 



12 A 



4.4 A 
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Perspective View of Fe Chain 




Electron Density Spin Excess 



FIG. 7: Ground state structure of Ci2oFeio. The axial posi- 
tions of the Fe atoms alternate between ±0.21 A. Also shown 
are the electron density cross-section (center) and spin excess 
isosurface (far right). 



Atomic Metal Necklaces within a C120 Torus 

The ground state was obtained for Mi loops encap- 
sulated within the C120 torus with the resulting struc- 
tures shown in Figures [7| [8] and [9j The associated prop- 
erties are summarized in Tables 1-3. Consistent with 
the ground state morphology of free standing metallic 
chains, the encapsulated necklace of Fe atoms forms a 
zigzag structure which is highlighted in the bottom of 
the figure. The minimum Fe-C bond length was found 
to be between that of the (3,3) and (4,4) SWNTs. The 
toroidal composite has a HOMO/LUMO bandgap of less 
than 0.1 eV and has a magnetic moment of 2.28/i# per 
atom-approximately the same as that associated with 
bcc Fe. The spatial distribution of this magnetic mo- 
ment is shown in Figure [7] through an isosurface plot of 
the spin excess. At zero temperature, the Fe atoms have 
aligned spin axes and so the composite is ferromagnetic. 
In contrast, neither of the straight nanotubes with Fe 
showed any intrinsic moment because of stronger Fe-C 
bonding for the tubes. 

Encapsulation of an Auio necklace was found to be 
energetically unfavorable, while the Cuio encapsulation 
is an energetically neutral process. As indicated in Tables 
2 and 3, these results are consistent with those associated 
with CNT encapsulations. 

The DFT investigation indicates that the C120 torus 
with all three types of metal necklaces is structurally 
stable at zero temperature. The TCN has a strong affin- 
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143° 

Perspective View of Au Chain 



FIG. 8: Ground state structure of C120AU10. The axial po- 
sitions of the Au atoms alternate between ±0.076 A Also 
shown are the electron density cross-section (center) and spin 
excess isosurface (far right). 




Perspective View of Cu Chain 



Structure 


M-M 


M-C 


M-M 


Gap 


Mom. 


E B 




A 


A min 


Angle 


eV 


I^b/M 


eV/M 


Au fee 


2.88 


- 


- 


0.00 


0.00 


3.05 


Au chain 


2.83 




57° 


0.00 


0.00 


2.07 


(3,3)-Au 


2.46 


2.34 


180° 


0.00 


0.00 


-6.80 


(4,4)-Au 


2.52 


2.56 


156° 


0.00 


0.00 


-0.12 


C120-AU10 


2.43 


2.00 


143° 


<0.1 


0.00 


-5.54 



TABLE II: Summary of properties for Au structures. The 
bandgap is HOMO-LUMO for the toroidal molecule. 



Structure 


M-M 


M-C 


M-M 


Gap 


Mom. 


Eb 




A 


Amin 


Angle 


eV 


Mb/M 


eV/M 


Cu fee 


2.56 






0.00 


0.00 


3.62 


Cu chain 


2.42 




59° 


0.00 


0.00 


2.17 


(3,3)-Cu 


2.46 


2.24 


177° 


0.00 


0.00 


-0.83 


(4,4)-Cu 


2.47 


2.69 


168° 


0.00 


0.00 


1.09 


C120-CU10 


2.39 


1.92 


144° 


<0.1 


0.00 


0.01 



TABLE III: Summary of properties for Cu structures. The 
bandgap is HOMO-LUMO for the toroidal molecule. 



FIG. 9: Ground state structure of C120CU10. The axial po- 
sitions of the Cu atoms alternate between ±0.025A Also 
shown are the electron density cross-section (center) and spin 
excess isosurface (far right). 



ity for Fe atoms (1.82 eV/Fe) and a neutral affinity for 
Cu atoms, but the placement of a Au necklace within 
the torus would require an energy input of 5.54 eV per 
Au atom. This precludes the study of Au wires within 
this setting. The (3, 3) and (4, 4) CNTs can be used 
to estimate the affinity of C120 toroids for metal atoms; 
the binding energy was found to be midway between the 
CNT values for each of the three metal chains. This 
allowed the prediction, in advance, that C120 would ex- 
hibit a neutral affinity for the Cu necklace. The pres- 
ence of the metal atoms causes the TCN composite to 
have a HOMO/LUMO bandgap of less than 0.1 eV in all 
three cases. When encapsulated with an Fe necklace, the 
torus exhibits ferromagnetic properties close to that of 



Structure 


M-M 


M-C 


M-M 


Gap 


Mom. 


E B 




A 


Amin 


Angle 


eV 


/xb/M 


eV/M 


Fe bec 


2.49 






0.00 


2.22 


4.89 


Fe chain 


2.35 




59° 


0.00 


3.25 


2.87 


(3,3) -Fe 


2.59 


1.94 


144° 


0.48 


0.00 


1.75 


(4,4)-Fe 


2.77 


2.22 


126° 


0.00 


0.00 


1.95 


Ci2o-Feio 


2.49 


2.04 


139° 


<0.1 


2.28 


1.82 



TABLE I: Summary of properties for Fe structures. The 
bandgap is HOMO-LUMO for the toroidal molecule. 



bec iron, and it will be useful to determine to what de- 
gree, if any, this persists at finite temperatures. It is not 
clear what sort of chain structure and composite proper- 
ties would result from the encapsulation of chains with 
symmetry at odds with that of the host torus. 

SUMMARY 

The azimuthal symmetry of encapsulated atomic 
metallic necklaces within toroidal carbon sheaths of- 
fers the prospect of synthesizing composites with novel 
electromagnetic properties that are not associated with 
straight, metal-encapsulated CNTs. As a first step in 
this direction, a particularly simple encapsulation was 
analyzed using density functional theory to determine 
its ground state geometry, structural stability, binding 
energy, electrical conductivity and magnetic moment. 
Of the four known classes of TCNs, a Type 1 C120 
molecule was chosen because it results in a computation- 
ally tractable composite, and only ten metal atoms were 
considered in order to take full computational advantage 
of the D 5 d symmetry of the torus. The TCN exhibits a 
positive, neutral and negative affinities for Fe, Cu and 
Au atoms, respectively, and metal atom encapsulation 
resulted in linearly stable structures in all three cases. In 
contrast to Fe chains within (3, 3) and (4, 4) CNT's, Fe 
atoms within the Ci2o-Feio molecule exhibit a net mag- 
netic moment comparable to bec Fe. 

All toroidal nanotubes should all be capable of serv- 
ing as sheathing hosts for atomic metal necklaces with 
novel electromagnetic properties. The possibility of en- 
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capsulating metal chains within the Ceo-coalesced tori of 
Types 3 and 4 is of particular interest because these are 
structures that have already been synthesized [6 ,[7 . For 
Type 3 TCNs, a reasonable strategy might be to encap- 
sulate Fe atoms within the Ceo molecules prior to coa- 
lescing them into a toroid. Sufficiently long CNTs with 
embedded metal atoms could also be deformed to syn- 
thesize Type 2 TCN composites. Both approaches would 
admit the construction of alloyed necklaces as well — the 
one-dimensional, atomic analog of powder metallurgy. 
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